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Abstract: We study the atom-light interaction in the fully quantum regime, with the 
focus on off-resonant light scattering into a cavity from ultracold atoms trapped in an 
optical lattice. The detection of photons allows the quantum nondemolition (QND) 
measurement of quantum correlations of the atomic ensemble, distinguishing between 
different quantum states. We analyse the entanglement between light and matter and show 
how it can be exploited for realising multimode macroscopic quantum superpositions, such 
as Schrodinger cat states, for both bosons and fermions. We provide examples utilising 
different measurement schemes and study their robustness to decoherence. Finally, we 
address the regime where the optical lattice potential is a quantum dynamical variable and is 
modified by the atomic state, leading to novel quantum phases and significantly altering the 
phase diagram of the atomic system. 
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1. Introduction 

Ultracold atomic systems in optical lattices form a powerful tool for studying and simulating the 
behaviour of quantum degenerate matter. Substantial interest in these systems has arisen partly as 
a result of their interdiseiplinary nature, as they provide exeiting possibilities to realise “toy-model” 
Hamiltonians, erueially, with highly tunable parameters, that ean deseribe a wide range of diverse 
physieal systems and proeesses from many different fields. Sueh phenomena inelude quantum 
phase transitions, strongly eorrelated materials, high-energy physies [1] and biologieal systems [2]. 
Additionally, they provide a strong potential avenue for performing quantum information proeessing 
protoeols. Light plays a fundamental role in these experiments, as it allows for eooling, trapping and 
manipulation of the atomie ensemble. Nevertheless, the light field is usually deseribed elassieally, while 
its quantum properties are left negleeted. However, elevating the treatment of light to the quantum 
regime leads to many interesting phenomena not present with elassieal light; pioneering works in the 
field have already shown that these phenomena inelude non-destruetive quantum measurement [3-8], 
self-organisation and other novel phase behaviour [9-15], as well as the possibility to engineer quantum 
states and dynamies, partieularly through the measurement baekaetion effeet [16-23]. 

In this artiele, we exploit the opportunities offered by the quantum nature of light for use in probing 
quantum eorrelations and realising maeroseopie quantum superposition states. To this end, we build upon 
earlier work on bosons [24], expand studies of quantum light as a quantum non-demolition measurement 
deviee and as a method of state preparation to the fermionie ease and make use of some novel possibilities 
that arise as a result of their spinful nature. We present a homodyne deteetion seheme, whieh may be used 
to engineer states, ineluding Sehrodinger eat states. Moreover, we also study some of the properties of 
the light-matter entanglement that forms a key ingredient in making sueh effeets and protoeols possible. 
Further, we show how the quantum eoupling between light and matter modifies the atomie trapping 
potential, leading to novel quantum phases and allowing for the implementation of tunable long-range 
interaetions not aehievable with 5-wave Feshbaeh resonanees. 

In analogy (yet eontrastingly) to eavity QED experiments [25], where the quantum light states in 
a eavity were probed and projeeted by measuring atoms, here, the roles are reversed, and we probe 
properties of matter fields with light. One sueh advantage this approaeh yields is in the ease of sealability; 
the number of atoms ean be tuned from a few to thousands by ehanging the number of illuminated sites, 
while sealing eavities is a ehallenge [25]. Thus, this work, alongside reeent experiments involving 
Bose-Einstein eondensates (BECs) in optieal eavities [26-28], where sueh predietions ean be tested, 
helps to bridge the gap between quantum opties and quantum gases, approaehing the fully-quantum 
regime of many-body light-matter interaetions [24,29] . 

We foeus on a system of ultraeold atoms trapped in an optieal lattiee with M sites and eoupled to 
quantum light. Speeifioally, the atoms are illuminated with an off-resonant eoherent probe beam with 
amplitude ao and frequeney cUp; the seattered light at a partieular angle ean be seleeted and enhaneed 
by an optieal eavity with deeay rate k and mode frequeney tUc (Eigure 1). The Hamiltonian of the full 
light-matter system is "H = deseribes the atomie dynamies in an optieal lattiee [1]. 

Eor bosons, this is the Bose-Hubbard Hamiltonian: 
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'^b = -to + h.c.) + Y 

(*.i> * 

and for fermions, the Hubbard model: 

■^Z = “^0 + ^-c.) + (2) 

o-e{t4} {id) i 

where bi and /j ^ are the bosonie and fermionic annihilation operators at site i, respectively, to is 
the nearest-neighbour hopping amplitude, U is the on-site atom-atom interaction energy, the bosonic 
(fermionic) number operator hi = = fi^fm) and the sum {i, j) is taken over neighbouring pairs 

of sites. The light field is described by = ha)ca\ai, where cii is the cavity mode annihilation operator. 
In the bosonic case, the light-matter interaction is given by [24]: 

= hQiid\aiFii + hQio ^agfiiT’/g -f aoojT’ig^ (3) 

where the atoms couple to the cavity mode via the coefficient ^lim = gigm/^a, gi is the light-matter 
coupling constant for mode /, is the detuning between light and the atomic resonance, Fim = Dim + 
Bim and: 

M M 

^ J'Tnu Bi,. = I] JirUh (4) 

{id) 


probe 



Figure 1. Setup. Atoms in an optical lattice are illuminated by a probe beam (red), and the 
scattered light (blue) is selected and enhanced by an optical cavity. 


The coefficients characterise the coupling between light and matter at each lattice site and are 
defined as the convolution between the light mode functions Mz(r) and the localised atomic Wannier 
functions w(r), i.e.. 


jlm _ 


w{y — ri)uKr)um{r)w{r — rj) dr 


( 5 ) 
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Note that the operator Dim deseribes the seattering from on-site density, while I3im that from small 
intersite densities [8]. For well-loealised atoms, JD ul{ri)uo{rj)Sij, and the eontribution to the light 
seattering due to Bim ean often be negleeted. 

In the fermionie ease, we introduee the light polarisation as an additional degree of freedom for the 
light field, whieh allows for probing of different spin speeies. In partieular, we use linearly-polarised 
light aix and aiy, whieh eouple differently to the two spin densities and beeause of seleetion 
rules that eonstrain the allowed transitions between different hyperfine states of the atoms [30-35]. 
In eontrast to the (spinless) bosonie ease, this allows us to probe different linear eombinations via 
the operators Dim,x = E* Jlr9i or Dim,y = J-f m* where p* = hi^ + hn (density) and 
= hi^ — Uii (magnetisation). In the fermionie ease, we also eonsider attraetive partiele-partiele 
interaetions. Additionally, due to the differing partiele statisties and Pauli exelusion, fermions exhibit 
different dynamies and phase behaviour to their boson eounterparts. 


2. Quantum Nondemolition Probing of Many-Body States 


The atom-light eoupling in the Hamiltonian introduees eorrelations between light and matter. 
As a eonsequenee, it is possible to eharaeterise the quantum state of the atoms by performing 
measurements of the light observables. Speeifieally, probing the light amplitude oi or photon number 
a\ai deseribes a quantum nondemolition (QND) measurement of the observables related to Dim [36]. 
We foeus on the regime where the atomie dynamies {i.e., tunnelling) is mueh slower than the light 
seattering: the tunnelling amplitude to and the interaetion U determine the quantum properties of the 
atomie state, but sueh properties do not ehange during the measurement time. The stationary solution of 
the Heisenberg equations for the seattered light operator oi (negleeting the small dispersive frequeney 
shift DiiDii) is given by: 


Oi 


HiotXo 

Ap + iK 


Dio 


CD 


( 6 ) 


where Ap = Wp — cUc is the probe-eavity detuning. Analogously to elassieal physies, the light amplitude 
(oi) depends on the mean density at eaeh lattiee site and deseribes the seattering of the probe beam ao 
from a periodie grating of atoms. Its angular distribution is determined by the position of the atoms: (cii) 
reaehes maximum value if they seatter light with the same phase, while it vanishes if the eontributions 
from eaeh lattiee site interfere destruetively. Moreover, quantum density-density eorrelations of the 
atomie state are imprinted on the expeetation value of the photon number operator: 


(ala,) = ichd'c) = (?> 

Higher moments of di probe higher-order eorrelation funetions (e.g., the fluetuations of the photon 
number depend on the four-point density expeetation values). In order to emphasise the differenee 
between elassieal and quantum seattering eontributions, we foeus on two-point eorrelation funetions 
and define: 


R = {£>'£)) - |(D)|= = ^ ((«,«,) - (fii)(a,)) 


( 8 ) 



Atoms 2015, 3 


396 



Figure 2. Classical diffraction patterns (normalised to iV^) and quantum additions 
(normalised to N) for N fermions in a one-dimensional optical lattice (45 sites, half filling) 
in the non-interacting regime (first row) and for U/to = 10 (second row). The lattice extends 
into the vertical direction, and the probe is along the horizontal, i.e., the angle marked with 
0° indicates the light intensity in the forward direction. The classical diffraction patterns 
(a,d) do not depend on the interaction while the quantum additions Ry (b,e) and (c,f) 
distinguish between different ground states. Note that the scale in (f) is different from (a-e): 
the fluctuations in the atom number increase as the attraction between the atoms favours 
doubly-occupied sites. Red circles indicate the orientation of the lattice. 

This quantity depends on the scattering angle via the coefficients Ju and is directly related to the 
fluctuations of the atom number. Therefore, quantum light scattering distinguishes between different 
quantum states, such as the superfluid and Mott insulator states [24]. As mentioned in the previous 
section, we extend the bosonic model [24,37,38] to also describe spin-1/2 fermionic atoms and define 
the quantum additions R^ and Ry, related to the fluctuations in density and magnetisation, respectively. 
If the mode functions of the probe and scattered light are travelling waves with wave-vectors kin and kout, 
respectively, the coefficients Jjj are proportional to exp [/(kin — kout) ■ r^], and the quantum addition 
R is proportional to a structure factor, i.e., the Fourier transform of the density-density correlations. 
Figure 2 compares the classical diffraction pattern and quantum additions Rx and Ry for fermions in 
a one-dimensional optical lattice at half filling. The scattering patterns were calculated for the ground 
state, obtained via imaginary time evolution using the TNTlibrary [39]. In particular, we describe the 
system using the Hubbard model and focus on the ground state of the system in two different regimes: 
non-interacting (U/to = 0) and strongly-attractive interactions (U/to = 10). Note that in both cases, 
the local magnetisation of the system is zero ((rhi) = 0), and as a consequence, the classical diffraction 
pattern for linear-?/ polarised light vanishes, i.e., {aiy) = 0. Nevertheless, the quantum addition Ry is 
non-zero and depends on the quantum state of the system. In particular, we find that Ry decreases with 
increasing values of U/to, as the attraction between the atoms favours doubly-occupied sites, i.e., the 
formation of pairs of fermions with opposite spin, which suppress the fluctuations in the magnetisation. 
Furthermore, the classical diffraction pattern due to the atomic density does not depend on U/to, as the 
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local density is independent of the interaetion for a translationally-invariant optieal lattiee ((pj) = 1). 
However, the presenee of doubly-oeeupied sites inereases the fluetuations in the atom number, leading to 
a stronger signal for R^. In general, the seheme we propose allows direet probing of the phase diagram 
for the superfluid-Mott insulator-Bose glass phases in the one-dimensional Bose-Hubbard model [8] 
and the deteetion of different quantum states of atomie systems sueh as dimers [14], density waves and 
magnetie order. 

3. Measurement Backaction and State Preparation 

Due to the entanglement between light and atoms, the deteetion of a photon affeets the dynamies 
of the atomie state. This is one of the most important manifestations of quantum meehanies, the 
so-ealled measurement baek aetion. In order to study this effeet, we foeus on the dynamies of the 
quantum state eonditioned to the measurement proeess in a single experimental realisation and deseribe 
it using the quantum trajeetory formalism. We eonsider fermions in a one-dimensional optieal lattiee, 
inside an optieal eavity with deeay Ky, sueh that only photons linearly polarised along the y axis are 
allowed to eseape it, and thus, the measurement seheme probes the magnetisation of the atomie sample. 
We illuminate K lattiee sites at equal intensity with a eoherent beam, and we deteet the ^/-polarised 
seattered light in the diffraetion maximum, effeetively probing the magnetisation of the illumined sites 
as aiy = C{Nk^ — Nri) = CMr, where Nr and Mr are the total atomie oeeupation number and 
magnetisation of all illuminated sites. If the i-th photon is deteeted at time ti, the state of the system 
ehanges instantaneously: the eavity annihilation operator aiy is applied to it as |'hc(fi)) —^ oiy l^c(fi)) 
(and the state is subsequently normalised). Moreover, the evolution of the system between two 
eonseeutive photo eounts is determined by the non-Hermitian Hamiltonian: 

Heff = -f — zKyalyOiy (9) 

The photo eount events and the evolution due to have opposite eontributions to the expeetation 
value of Mr, as the former tends to inerease the magnetisation of the illuminated area, while the latter 
tends to deerease it. The full dynamies of the system is therefore eonditioned by the photodeteetions 

oeeurring at times ti,t2 ., whieh are determined stoehastieally. The initial wavefunetion of the 

light-matter system ean be written: 

l«'(0)> = lqt) |kO Kk„(0)) Kki(0)> (10) 

k,q 

where |(Xqkp(0)) are the eoherent states of the light field with polarisation p (in the eireular polarisation 
basis) and |qcr) are the Foek states for fermions of speeies cr. Sinee we are negleeting atomie tunnelling, 
the Hamiltonian H®® does not mix the different atomie Foek states, and it is possible to write an analytieal 
expression for the eonditional evolution of |'kc(f) )• Speoifieally, the eoherent states |aqkp(0)) aequire 
a phase faetor exp[4>qk(f)], whieh depends only on the atomie eonfiguration qk [18,19,24,40]: this 
phase reaehes a steady state value after f 3> 1/ for eaeh quantum trajeetory, allowing us to write the 
eonditional wavefunetion of the system after m photodeteetion events as: 

) cx (oCkqi? - (Xqki)"* |q t) |k |) |aqkfl(f)) \ct^kL{t)) 

k,q 


( 11 ) 
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Due to the steady state in all of the light amplitudes (Xqkp, this expression does not depend on the 
speeifie deteetion times U. Interestingly, the state of the system is not faetorisable into a produet of matter 
and light states as the two remain entangled. Moreover, the eavity introduees an effeetive interaetion 
between the two spin speeies, sinee different atomie eonfigurations with different spins are eoupled 
as a result of the photodeteetion. Foeussing on the properties of the atomie system, we eompute the 
probability distribution of having fermions with spin cr in the illuminated area of the optieal lattiee at 
time t. This ean be obtained summing the absolute values of the eoeffieients of all of the eonfigurations 
q'k' with the same and introdueing the initial probability distribution P^^\za) = Xlq'k' I'^q'k'P- 
Henee, the eonditional probability distribution is: 

Pc{t,z^,zi) = - ^J^™exp (-r|zt - P^^\z^)P^^\zi) (12) 

where AT is a normalisation eonstant and r = 2\C\‘^Kyt. Rewriting this expression in terms of the 
magnetisation of the system, we have: 

Pc{t,MK) = ^\MK\^^exp pW(^ir) (13) 

The relation between m and t follows a stoehastie proeess, whieh defines a single quantum trajeetory. 
The measurement proeess strongly modifies (Mk)' when the first photon is deteeted, this probability 
distribution beeomes bimodal, and all of the eonfigurations with = 0 are forbidden (Figure 3). 
Moreover, as time progresses, P(f, Mk) shrinks to two narrow peaks around the values Mi 2 = ±\/m/T 
with deereasing width (note that the stoehastie nature of this expression arises through m). Thus, the 
final state of the system has a well-defined absolute magnetisation as the measurement proeess projeets 
the atomie state to a superposition of two states with magnetisations Mi and —Mi: a Sehrodinger oat 
state. This multioomponent struoture emerges as a oonsequenoe of the degeneraey of the photon number 
operator d\ydiy and ean be exploited to realise non-trivial quantum superpositions by deteoting photons 
at different angles [24] or by engineering the eoeffieients Ja using different light mode funetions [22,23]. 
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Figure 3. Probability distribution P{z^, z^) without measurement (a) and with measurement 
(c) for a fermionie system. The insets show the probability distribution of the magnetization 
P{Mk) for the initial state (b) and after photons are deteeted (e). Measurement baekaetion 
ereates a Sehrodinger eat state. 
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4. Light-Matter Entanglement 

The light-matter entanglement is a key feature of the system; indeed, it is this entanglement that is 
responsible for many of the effects observed when one goes beyond semiclassical treatments. It is thus 
interesting to consider the properties of this entanglement and how it is affected by the measurement 
backaction. As we deal with pure states, we can characterise the (bipartite) entanglement properties with 
the entanglement entropy [41], defined as the Von Neumann entropy of the reduced density matrix of 
one of the subsystems: 

E{\^)ab) = S{9a) = -Tr(pAlogp^) (14) 

If the reduced density matrix is diagonal, this is equivalent to the Shannon entropy i?(p) = 
— Pi log Pi, where Pi are the diagonal elements of p. When the number of atoms is large, we can 
approximate the distribution of subsystem occupation numbers as continuous, such that we may invoke 
the central limit theorem and treat them as Gaussian functions [42]. For such functions, the Shannon 
entropy is given by: 

P(P) = ilog(27reo-2) (15) 

where cr^ is the variance of the distribution. 

We first consider the entanglement present in the steady state for bosons in a lossless cavity. In this 
case, when the lattice is illuminated by a coherent probe, the state of the system can be written |\l'(f —>■ 
cxd) ) = Cq |q) IcVq), and hence, the reduced state of the light as: 

PL ^ ^ I Cq I I OC-q ) ( CXql ( 16 ) 

q 

Thus, the entanglement entropy can be obtained directly from the Shannon entropy of the distribution 
-P(^) = Sq|aq=az I'^ql^’ where 2 ; represents all states q corresponding to the same cCq {i.e., z is the 
light measurement eigenvalue), and we have assumed that light states corresponding to different z are 
approximately orthogonal. We justify this by noting the overlap of two coherent states is given by 
I (o-i I(X 2 ) P = and, hence, is exponentially suppressed by the difference in the (X, which can be 

tuned by, for example, increasing the pump power. 

We investigate two paradigmatic examples where there is one incoming probe beam and one scattered 
beam. In the first, we illuminate a region of the lattice at the diffraction maximum (z = N^) and the 
second at the diffraction minimum (z = When the matter is initially in a superfiuid state, 

for the diffraction maximum, we have that 2 ; obeys a Poisson distribution P{z) = {NY/z\, where 

{N) is the expected number of atoms in the illuminated region, while for the diffraction minimum, it 
obeys a Skellam distribution P{z) = e~^^Yz{2{N)), where Y is a modified Bessel function of the first 
kind, and we have assumed that the expected number of atoms on odd and even sites is equal. Since both 
of these distributions have the same variance cr^ = (N), from Equation (15), we find that for both cases, 
the light-matter entanglement entropy is: 

E = ^log{27re{N)) (17) 

When cavity loss is taken into account, the Cq and, hence, P{z) evolve, according to P{z,t) = 
|p^m| 2 g- 2 |a 2 pKtp(^ g)/^ [24]. TMs squcczcs the initial distribution, and hence, the light-matter 
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entanglement will deerease over time and eventually vanish, as the distribution of 2 ; eonverges to a single 
value; henee, the system enters a produet of the light and matter states. 

However, a peeuliarity arises in the diffraetion minimum ease if we eonsider only the leaked photon 
eount rate to have been deteeted and the phase not measured. In doing so, only the magnitude of is 
eonstrained, and there is no distinetion between 2 ; and — 2 ;. Thenee, at long times, the state of the system 
beeomes a eat state [24] : 

|^(r,m ^ oo))^(|^) |(X^) + (-1)”^ \-z) |(X_^)) (18) 

possessing a light-matter entanglement entropy of one independent of the number of atoms, 2 ;, and the 
light field strength. 


5. Homodyne Detection 


The measurement seheme introdueed in the first seetion of this artiele is extremely flexible and ean 
be easily extended. For example, we eonsider a setup that ineludes a loeal oseillator eavity and another 
eavity eontaining the atomie system, as illustrated in Figure 4, where k and y are eavity deeay eonstants 
for the atomie system and loeal oseillator, respeetively, and (3 is the loeal oseillator eoherent amplitude. 
We eombine the two light modes with a beam splitter with refleetivity R, and in order for all of the 
atomie eavity photons to reaeh the deteetor in our system, we let i? —)■ 0. However, this also requires 
that the loeal oseillator amplitude, |3 —)■ 00 and y —)■ 0, sueh that the loeal oseillator flux at the 
deteetor, 3F = i?^y|(3p, remains eonstant [43]. Assuming that the atomie dynamies is mueh slower 
than the light seattering, it ean be shown that the evolution of the light-matter system is deseribed by the 
non-Hermitian Hamiltonian: 

Heff = R} A- - ihKa\ar - (19) 

where 0 is the loeal oseillator phase at the deteetor. In this ease, the jump operator that is applied to the 
quantum state at eaeh photodeteetion is [43]: 

c = (20) 


The state of the system at time t after m photodeteetions is: 

\'^c{m,t)) oc c°_^ [ 1 -f -t-c°_ [l-f e*^'^ 2 ;_]”'e"*‘^*| 2 ;_)|a^_) ( 21 ) 


where 93 is a eonstant, whieh ean be derived from the Hamiltonian, 


z± 


c± 




sin^(60) 


cos(60) 


( 22 ) 

(23) 


where 60 = — 0 is the phase differenee between ao and the loeal oseillator, C = 

and 2 ; are the unique eigenvalues of the D operator, and as sueh, they are a suitable label of the 

measurement eigenstates. 
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Figure 4. The homodyne detection scheme as analysed in [43]. The pumped cavity contains 
an atomic gas in an optical lattice potential. The second cavity radiates a coherent local 
oscillator field ||3). 


The homodyne measurement scheme allows us the freedom to choose the phase difference between 
the local oscillator and the photons emitted from the atomic cavity system, bcf) = (pc — 9, by adjusting 
the local oscillator phase. There are two interesting special cases depending on the value of this phase 
difference. Firstly, we can set the difference to an integer multiple of 27r, bp = 2nn. In this case, the 
final state is given by: 


\'^cim,t)) (X cly'^^\z+)\oc^^) +c°_(-l)'"e (24) 

with = ^yl/2K\C\'^{±^ym/t—'/T), which is a cat state, albeit still entangled with the light. The light 
and matter can then be disentangled, though, by switching off the probe and detecting all leaked photons, 
leaving the cavity in the vacuum mode [18]. However, each photo count flips the phase difference 
between the two components by tt, making this state particularly fragile with respect to decoherence: 
with each missed photo count, the state loses its purity and becomes a mixed state, losing any atomic 
entanglement it possessed. Another notable case is if the two phases are offset by 7r/2, bp = 7r/2 + 
27m. In this case, 2 ; = 2 ;^ = — 2 ;_ = ^/l/2K\C\‘^^/m/t — T, and the state evolution conditioned to the 
measurement leads to: 


(X Vt + iy/mjt — T 


jipt 


z)\oiz 


+ cL 


Vt — i\J mjt — T 




- 2 ; - a. 


(25) 


Importantly, mjt — T small, and each photo count imparts only a small change in phase between 
the two components. The state Equation (25) is therefore more robust than Equations (18) and (24), as 
it does not suffer much from the decoherence problem associated with a missed photo count, and this 
scheme may prove a more useful method to prepare cat states [18]. 
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6. Quantum Optical Lattices 

The quantum properties of light also induee a modifieation to the optical lattice potential experienced 
by the atoms. This leads to an effective quantum potential that changes the potential energy surface: the 
optical lattice becomes a quantum dynamical variable that depends on the atomic state. The effective 
Hamiltonian in the steady state of light is [12,14,44]: 

n^^ = n^ + g,sD^D (26) 

where ges = hAp\C\‘^. The cavity introduces an effective long-range interaction, which allows the 
system to support novel emergent orders, even in a single-mode cavity [14]. Therefore, the properties 
of the ground state of do not depend only on the local processes (tunnelling and on-site interaction) 
present in the usual (Bose) Hubbard Hamiltonian, but also on the detuning Ap via the £)-operator. This 
drastically changes the effective steady state of the system. It is instructive to analyse the simplest 
measurement scheme, where atoms scatter light homogeneously into a single-mode cavity. In this case, 
there is no spatial variation in the density couplings Ja = Jd, and thus, there is no density wave 
instability. Without quantum light and in the absence of disorder, the system supports superfluid and 
Mott-insulating phases [1]. However, the quantum potential will distort the phase diagram, and it is able 
to shift the superfluid-Mott insulator transition conditional on the light amplitude. 

Here, we show that light favours the superfluid state and stabilises it, even in deep optical potentials. 
Specifically, for a given value oHq/U corresponding to a Mott insulator state in the absence of quantum 
light, we find that the ground state of the system becomes a superfluid state with minimal atomic density 
fluctuations per site, i.e., only two Fock-state components, with integer fillings m and m + 1 per site. 
Therefore, the local atomic density p = (n*) is constrained between these two values: m < p < m -I- 1. 
Considering the case of a very deep optical lattice, where tunnelling can be neglected, the effective 
on-site energy can be written as: 

Est = U{g-l)g/2-Ugx-Up‘^/{2yD) (27) 

withy^) = U/{2K(Xd), ocd = g = x + I, x = {[i/U — p/yr?) and qis the chemical potential. 

Expressing Equation (26) in the Fock state basis, we can diagonalise it and self-consistently derive the 
on-site density p = (?).*), order parameter = (bi) and atomic fluctuations A(nj) that minimise the 
ground state energy for a homogeneous system [45,46]. The system is a superfluid provided that "0 7 ^ 0, 
and if not, the system is in the Mott insulator state. Specifically, we find that for miyo + 1 ) < \^/u < 
y))^ {m+l)+m, corresponding to the minimal fluctuation superfluid state, these parameters are given by: 


ijj = \/ {m + l){p — m){l — p + m) 

(28) 

1 

Q. 

(29) 

,.) = (p - m)(l - p -1- m) 

(30) 


while for y))^(m -f 1) -f m < \i./U < {m + l)(y))^ -f 1), one has ■^ = 0, p = m-|-l and A(fij) = 0. 
The presence of such a state can be easily seen in the mean occupation p, where in between the Mott 
insulator plateaus, it behaves linearly as a function of the chemical potential. The slope of this function 
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depends on the coupling and, in the large K limit, is The phase diagram in terms of the effective 
light strength oto shows that Mott insulator states are indeed strongly suppressed (see Figure 5), and 
in between the Mott insulator plateaus, these minimal fluctuation superfluids emerge. This behaviour 
is confirmed by exact diagonalisation simulations where, in the limit of large numbers of illuminated 
sites, the same pattern is evident. Moreover, we find that p(p,/[/) has discrete steps depending on the 
number of illuminated sites K and that they become smooth as iT —)■ oo. This computation can easily 
be implemented by noting that since we are focussing on the limit to/U = 0 and in absence of a trapping 
potential (besides the classical optical lattice), the number of basis states needed to find the ground state 
reduces significantly. Thus, one can easily achieve simulations with 100 or even 1000 sites, and at this 
point, the continuum limit is evident. We use the symmetries of the effective light-induced long-range 
interaction in a truncated Hilbert space where the maximum filling per site is five particles; as for the 
range of the on-site interactions and chemical potential considered, this is sufficient. 



Figure 5. Phase diagram of the bosonic system with quantum optical lattice for 
homogeneous scattering from the atoms in the limit 1^/11 = 0. (a) Superfluid order 
parameter: bright areas correspond to the superfluid state with minimal atomic fluctuations 
and dark areas to Mott-insulator regions, (b) Local density: Mott insulator plateaus surround 
minimal uncertainty superfluid states with consecutive Fock-state fillings, and white lines 
denote phase boundaries. 


7. Conclusions 

In summary, we have shown how to use light scattering from ultracold atoms to perform QND 
measurements on fermionic systems and demonstrated how the quantum addition to the classical 
diffraction pattern carries information about quantum correlations of the atomic state. Focussing on a 
single experimental run, we discussed the entanglement properties between light and matter and showed 
how these can be exploited for creating macroscopic quantum superpositions with fermionic and bosonic 
systems. We provided different schemes for the realisation of such states and suggested how to make 
them more robust to decoherence by using a homodyne measurement. Finally, we described the phase 











Atoms 2015, 3 


404 


diagram of a system of bosons in an optical cavity trapped in an optical lattice. The cavity mediates an 
effective long-range interaction that stabilises a superfluid state with minimal fluctuations and suppresses 
the Mott insulator states. 
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